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ABSTRACT 

General-purpose finite element computer codes that can model inelastic material 
behavior have been available for more than a decade. However, these codes 
have not been accurate enough for use in analyzing hot-section engine compo- 
nents. To correct this problem, General Electric developed a series of nine 
new stand-alone computer codes for NASA. Because of the large temperature 
excursions associated with hot-section engine components, these codes have 
been designed to accommodate broad variations in material behavior, including 
plasticity and creep. The capabilities of these computer codes are summarized 
here. 


*Work performed on-site at the Lewis Research Center for the Structural 
Mechanics Branch. 
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OBJECTIVE 


Under a two-year program at the General Electric Company, a series of three- 
dimensional inelastic structural analysis computer codes were developed and 
delivered to NASA. The objective of this program was to develop analytical 
methods capable of evaluating the cyclic time-dependent inelasticity that 
occurs in hot-section engine components. Because of the large temperature 
excursions associated with hot-section engine components, the techniques devel 
oped had to be able to accommodate broad variations in material behavior, 
including plasticity and creep (McKnight et al., 1986). To meet this objec- 
tive, a matrix consisting of three constitutive models and three element 
formulations was developed. A separate program for each combination of consti 
tutive model and element formulation was written, making a total of nine pro- 
grams. The source codes of the nine programs range in size from 7300 lines 
for the Bodner 20-node code to 19 000 lines for the Haisler-Allen nine-node 
code. The table below shows the length of each source code. All of the codes 
were given a stand-alone capability of performing cyclic nonlinear analysis. 


Constitutive 

models 

Element formulation 

20-Node 

8-Node 

9-Node 

Simple 

8300 

13,800 

17,900 

Classical 

9200 

16,300 

19,000 

Unified 

7300 

13,800 

17,600 


CD-86-31650 
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CONSTITUTIVE MODELS 


The three constitutive models are a simple model, a classical model, and a 
unified model. In an inelastic analysis the simple model uses a bilinear 
stress-strain curve to determine the plastic strain and a power law equation 
to obtain the creep strain. The second model is the classical model of Allen 
and Haisler (1981). The third model is the unified model of Bodner, et al. 
(1979). The attributes of the three constitutive models are shown below. All 
of the models were programmed for a linear variation of loads and temperatures, 
with the material properties being temperature dependent. 


Classical Unified 

• Uncoupled plasticity • Coupled plasticity 

and creep and creep 

• Plasticity 

Combined isotropic and Isotropic hardening 
kinematic hardening 

Piecewise linear No yield surface 

stress-strain curves 

Modified Prandtl-Reuss 
flow rule 

• Creep 

Steady state 
Isotropic hardening 
Prandtl-Reuss flow rule 
Nonisothermal 

CD-88-31663 


Simple 

• Uncoupled plasticity 
and creep 

• Plasticity 

Isotropic hardening 

Piecewise linear 
stress-strain curves 

Prandtl-Reuss flow rule 

• Creep 

Steady state 
Isotropic hardening 
Prandtl-Reuss flow rule 
Nonisothermal 
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ELEMENT FORMULATIONS 


The three element formulations available are an 8-node isoparametric shell ele- 
ment, a 9-node shell element, and a 20-node solid element. Both of the shell 
elements are obtained by degenerating three-dimensional isoparametric solid 
elements and then imposing the necessary kinematic assumptions in connection 
with a thin shell. The eight-node element uses serendipity shape functions 
for interpolation and Gaussian quadrature for numerical integration. The 
nine-node element uses Lagrange shape functions and Simpson's rule for numeri- 
cal integration. The 20— node solid element uses Gaussian quadrature for inte- 
gration. The attributes of these elements are listed below. 


8-Node shell 

• Five DOF 

3 displacements 
2 rotations 

• Serendipity shape 
functions 

• No rotational stiffness 
about normal to mid- 
surface. Deleted before 
assembly 

• Isotropic or 
orthotropic elastic 
properties 

• Surface, line, nodal, 
rotational, thermal, and 
gravity loads 

• Prescribed displacements 

• Gaussian quadrature 
used for numerical 
integration 


9-Node shell 

• Five DOF 

3 displacements 
2 rotations 

• Lagrange shape functions 

• Rotation about normal to 
midsurface treated as a 
prescribed displacement 

• Isotropic or 
orthotropic elastic 
properties 

• Surface, line, nodal, 
rotational, thermal, and 
gravity loads 

• Prescribed displacements 

• Simpson’s rule used for 
numerical integration 


20-Node solid 

• Three DOF 
(3 displacements) 


• Isotropic or 
orthotropic elastic 
properties 

• Surface, nodal, 
rotational, thermal, and 
acceleration loads 

• Prescribed displacements 

• Gaussian quadrature 
used for numerical 
integration 

CD-88-31664 
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INPUT FEATURES AND ANALYSIS TECHNIQUES 


To analyze structures with linear material behavior, the nine codes use a 
blocked-column-skyline, out-of-core equation solver. To analyze structures 
with nonlinear material behavior, the codes use an initial stress iterative 
scheme. To increase the convergence rate of the iterative scheme, Aitken's 
acceleration scheme was incorporated into the codes. 

The ability to model piecewise linear load histories was written into the 
codes. Since the inelastic strain rate can change dramatically during a linear 
load history, a dynamic time-incrementing procedure was included. The maximum 
inelastic strain increment, the maximum stress increment, and the maximum rate 
of change of the inelastic strain rate are the criteria that control the size 
of the time step. The minimum time step calculated from the three criteria is 
the value that is used. 

In dynamic analysis the eigenvectors and eigenvalues can be extracted by using 
either the determinant search technique or the subspace iteration method. 

These methods are only included with those finite element codes containing the 
eight-node shell element. The table below lists the features of the nine 
codes. 


Feature* 

r 



Model 



□ 


Simple 

Classical 

Unified 




Number of nodes 




8 

9 

20 

8 

9 

20 

8 

9 

20 

Free format data input 
Global coordinate system 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Cartesian 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Spherical 

X 

X 

na 

X 

X 

na 

X 

X 

na 

Cylindrical 

Local coordinate system 

X 

X 

na 

X 

X 

na 

X 

X 

na 

Cartesian 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Spherical 

X 

X 

na 

X 

X 

na 

X 

X 

na 

Cylindrical 

X 

X 

na 

X 

X 

na 

X 

X 

na 

Automatic generation of nodal coordinates 

X 

X 

na 

X 

X 

na 

X 

X 

na 

Automatic generation of element 
connectivities 

X 

X 

na 

X 

X 

na 

X 

X 

na 

Restart capability 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Dynamic allocation 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Biocked-column-skyline equation solver 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Initial stress iterative scheme 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Aitken's acceleration scheme 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Dynamic time incrementing 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Dynamic analysis 

X 

na 

na 

X 

na 

na 

X 

na 

na 

Material change option 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Element removal option 

X 

X 

na 

X 

X 

na 

X 

X 

na 

Midside node generation 

X 

X 

na 

X 

X 

na 

X 

X 

na 

Skewed coordinate system 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Orthotropic orientation definition 

X 

X 

X 

X 

X 

X 

X 

X 

X 


a x denotes "feature of the code"; na denotei 


"feature not yet available.” 


CD-88-31665 
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CYLINDRICAL SHELL ROOF 


The cylindrical shell roof shown below has frequently been used to test the 
behavior of shell elements. The eight-node element is compared with the four- 
node and eight-node elements of both MARC and NASTRAN. As can be seen from the 
graph the specialty eight-node finite element compares very favorably with the 
other elements. 


Z,W 



NORMALIZED 

DISPLACEMENT 



MARC (FOUR NODES) 

MARC (EIGHT NODES) 
NASTRAN (FOUR NODES) 
NASTRAN (EIGHT NODES) 
SPECIALTY FINITE ELEMENT 
(EIGHT NODE) 

ANALYTICAL SOLUTION 


ELEMENT ASPECT RATIO 


CD-B8-31666 
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